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We present a refined calculation of the quark-loop contribution to hadronic light-by-light scatter- 
ing that focuses upon the impact of the transverse components of the quark-photon vertex. These 
structures are compared and contrasted with those found within the extended NJL-models. We 
discuss similarities and differences between the two approaches and further clarify the important 
role of momentum dependent dressing functions. As expected we find that the transverse structures 
of the quark-photon vertex lead to a suppression of the quark-loop contribution to the anomalous 
magnetic moment of the muon. However, we find evidence that this suppression is overestimated 
within models with simple approximations for the quark-photon interaction. 



I. INTRODUCTION 

The anomalous magnetic moment of the muon is 
an observable that furnishes a precision test of the elec- 
tromagnetic (EM) interaction, the weak interaction and 
the strong interaction. The relative size of the different 
theoretical parts and comparison with the experimental 
results of the E821 experiment at Brookhaven [1, 2] are 
shown in table I. The dominant contribution of more 
than 99% is due to the EM interaction as described by 
quantum electrodynamics (QED), which has been eval- 
uated up to order a 5 in the fine structure constant [ ]. 
In addition, given the theoretical and experimental pre- 
cision available, the weak interaction yields significant 
contributions [4]. An interesting feature of the observ- 
able a M is its sensitivity to non-perturbative QCD correc- 
tions. The leading QCD contribution amounts to about 
0.006% which is about two orders of magnitude larger 
than the total theory uncertainty. Together, the leading 
and sub-leading QCD contributions dominate the error 
of the total Standard Model prediction as can be seen 
in the table. This is because perturbation theory is not 
applicable and the methods used lead to a substantial er- 
ror here. The two most relevant hadronic contributions 
are the hadronic vacuum polarisation (HVP) and the 
light-by-light (LBL) scattering contribution. The former 
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QED 


116 584 718.1 ( 0.2) 


99.99390% 


00.00098% 


weak 


153.2 ( 1.8) 


00.00013% 


00.07910% 


QCD LOHVP 


6 949.1 (58.2) 


00.00596% 


82.69628% 


QCD HOHVP 


-98.4 ( 1.0) 


00.00008% 


00.02441% 


QCD LBL 


105 (26 ) 


00.00009% 


16.50391% 


Standard Model 


116 591827.0(64 ) 


100% 


100% 


Experiment 


116 592 089 (63 ) 






Exp-Theo 


262 (89 ) 







TABLE I. Standard Model contributions to the muon g — 2. 



FIG. 1. The light-by-light scattering contribution to the muon 
g — 2. The main ingredient is the hadronic photon four-point 
function ILpvap to be discussed below. 



is related to experimentally available e + e~ — > hadrons 
data via dispersion relations (see [5]) and as such the 
error can be systematically reduced. This quantity is 
furthermore the subject of several recent lattice stud- 
ies [ ]. Furthermore, we have applied the method of 
Dyson-Schwinger equations (DSEs) to this quantity and 
were able to reproduce the dispersion analysis result on 
the ten percent level [9]. The result for HVP quoted in 
Table I is taken from Ref. [10]. Whilst presently the un- 
certainty is dominated by HVP, the LBL contributions 
(shown in Fig. 1) are potentially more problematic in the 
long run since it is extremely hard to determine these 
in a model-independent way. The LBL contributions 
have been investigated from the view point of low-energy 
effective models such as the Extended Nambu-Jona- 
Lasinio (ENJL) model [11], the Hidden Local Symme- 
try (HLS) model [12], vector meson dominance (VMD) 
approaches [13, 14], the non-local chiral quark model 
(NLxQM) [15, 16], the chiral constituent quark model 
(xGQM) [17], in holographic models [18] and Dyson- 
Schwinger Equations [19, 20]. The lattice calculations 
of LBL are still at an exploratory stage [2 1 ] . The LBL- 
contribution quoted in Table I is taken from Ref. [22]. 
There different groups, pursuing the strategy of hadronic 
models, agreed on this number. 

A future experiment, to be conducted at Fermilab, will 
measure the anomalous magnetic moment of the muon a M 
to a precision of 0.14ppm [23]. It is thus mandatory to 
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work towards getting the LBL contribution under suffi- 
cient control. For this undertaking we require mature 
non-perturbative methods that are well- rooted in QCD, 
such as DSE's and lattice QCD. We believe that a promis- 
ing way for the future is to combine these methods in a 
complementary fashion. 

In previous works we have determined important parts 
of the hadronic LBL contributions such as pseudoscalar 
meson exchange and non-transverse contributions to the 
quark- loop part of LBL [19, 20]. What has been left out 
so far are the transverse structures of the quark-photon 
coupling due to numerical complexity. From vector- 
meson dominance models, however, these contributions 
are believed to be sizable and negative and therefore lead 
to a substantial overall reduction of the LBL contribution 
to the muons anomalous magnetic moment. A complete 
approach towards LBL therefore has to include these con- 
tributions explicitly. In this work we provide a further 
step into this direction. The main focus of this article is, 
however, on the comparison between our approach and 
models like the ENJL or chiral quark model. We will ar- 
gue that such models provide simple tools to give impor- 
tant qualitative insights into the significance of different 
contributions to LBL. However, in order to provide more 
precise quantitative results, more elaborate approaches 
that take into account the full momentum dependence of 
dressing functions are mandatory. 

This paper is organized as follows. First we briefly in- 
troduce the quantity under question and give the Dyson- 
Schwinger Equations (DSE) and the particular trunca- 
tion used in this work, together with some notation, in 
section II. In section III we introduce the hadronic four- 
point function that lies at the heart of the LBL contri- 
bution and define it in terms of Green's functions. The 
main body of this work, the quark loop contribution to 
LBL, is discussed at length in section IV. Here we make 
comparison between the DSE and the ENJL approach in 
order to highlight similarities and differences. This is fol- 
lowed by a detailed discussion of our numerical results in 
section V, where the main focus is on the influence of mo- 
mentum dependent dressing functions and the different 
structures in the self-consistent quark photon vertex that 
are a vital part of our calculation. Details concerning the 
quark-photon vertex and the derivation of the hadronic 
four-point tensor in the present truncation are given in 
the appendix. 



the muon mass-shell this vertex can be decomposed as 




u(pt) 



2m,, 



nip), (1) 



where p and p' are the muon momenta, q is the photon 
momentum and <r a p = | [7^,7/3]. The anomalous mag- 
netic moment is defined as 
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«m = V- = *a(0), 



(2) 



which is obtained from Eq. (1) in the limit of vanish- 
ing photon momentum, q 2 . Here we will use the tech- 
nique advocated in Ref. [24] which simplifies the numer- 
ics by ensuring that all integrals are explicitly finite, see 
Ref. [20] for details. 

In the following we introduce the Dyson-Schwingcr 
equations for the quark propagator and the quark-photon 
vertex, together with the truncation scheme used in our 
calculations. The dressed quark propagator is given by 



S- 1 (p) = Zj 1 (p 2 ){-ijp' + M{p 2 )) 



(3) 



where Zf(p ) is the quark wave-function renormalisation 
and M(p 2 ) is the quark mass function. These scalar func- 
tions are obtained as solutions to the quark DSE, 



S- 1 (p) = Z 2 S 



g 2 Z lf ^ J dk^S(k)T"(k,p)D^(q) 



where dk = d 4 k/ '(27r) 4 and q = k 



(4) 

p is the gluon's mo- 



mentum. The bare inverse quark propagator is Sq (p) = 
—iff + rrio- This bare mass is related to the renormalized 
one by ZtfriQ — ZiZ m m q , with Z2 and Z m the wave- 
function and quark-mass renormalisation constants. The 
renormalisation constant for the quark-gluon vertex is 
Zif. To solve Eq. (4) we need the gluon propagator 
D^v{q) and the quark-gluon vertex T v {k,p). 

The quark-gluon interaction that appears in the DSE 
for the quark reads: 



(5) 



In Landau gauge, = T llv (q)Z{q 2 ) / q 2 where the 

transverse projector is T^^q) = 5^ v — q^q v /q 2 . The 
quark-gluon vertex T v (k,p) can be decomposed into 
twelve Dirac covariants. However, we will employ the 
rainbow-ladder (RL) truncation, which requires that we 
replace the complicated structure of the quark-gluon ver- 
tex with just its 7^ component. Hence Eq. (5) becomes 



II. BASICS 



(6) 



To obtain the LBL contribution to the muon anoma- 
lous magnetic moment a^ BL one must consider its contri- 
bution to the muon-photon vertex shown in Fig. 1. On 



where A(g 2 ) is the non-perturbative dressing of the 7„ 
part of the quark-gluon vertex, restricted to depend only 
on the exchanged gluon momentum. Combining all scalar 
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dressings into one effective running coupling, a e g(q 2 ) we 
have 



■lu, (7) 



where a e ff is a renormalization group invariant. The fac- 
tor Z\ ensures multiplicative renormalizability. 

In contemporary Dyson-Schwinger studies one employs 
the Bethe-Salpeter equations to study bound states of 
two particles, with the interaction described by a two- 
body kernel. To provide a realistic description of pseu- 
doscalar mesons one requires that the dynamical break- 
ing of chiral symmetry is encoded into the truncation. 
The symmetry-preserving two-body kernel correspond- 
ing to this yields the 'ladder' part of RL. For simplicity 
we quote this in terms of a c g, 



K rs ,tu(q) = 47rZ; 



a e fi(q 2 



(8) 



Note that it is possible to employ a beyond-RL truncation 
here, see Refs. [25-29]. In practice such an extension 
complicates the numerics considerably and is therefore 
not yet viable in the context of g — 2. 

We will employ an effective interaction called the 
Maris- Tandy (MT) model which has much phenomeno- 
logical success for pseudoscalar and vector meson masses, 
decay constants and form factors [. -34] . Success in the 
meson sector has led to its widespread use in the calcula- 
tion of baryon properties [35-38] . This effective running 
coupling is given by 



D 



x 2 e 



2 / 2 

q /uj . 



y 



27r 7m (l-e-") 
' log [e 2 - 1 + (1 + z) 2 ] ' 
q 2 /A 2 , z = q 2 /A 2 QCD 



(9) 



and features a Gaussian distribution in the infrared that 
provides dynamical chiral symmetry breaking. It is char- 

1/3 

acterized by an energy scale (wD) = 0.72 GeV, fixed to 
give the pion decay constant, and we choose u = 0.4 GeV. 
The second part reproduces the one-loop running cou- 
pling at large perturbative momenta. It includes the 
anomalous dimension j m = 12/(lliVc — 2iV/) of the 
quark propagator, and we use 7 m = 12/25, Aqcd — 
0.234 GeV and A f = 1 GeV. Note that we also employ a 
Pauli-Villars like regulator with a mass scale of 316 GeV. 
We focus here on the two lightest quarks whose mass at 
fi = 19 GeV is 3.7 MeV. 

The equation for the quark-photon vertex can be writ- 
ten in the form of an inhomogeneous Bethe-Salpeter 
equation in rainbow-ladder truncation, dependent upon 
the same two-body ladder kernel Eq. (8) 



[T M (P,fc)]„ = Z i7m 



(10) 



~ Z 2g / dq[S(q + )T tl (P,q)S(q-)} ut Ktu,rs(k - q) , 

where dq = <i 4 g/(27r) 4 , P is the outgoing photon mo- 
mentum, q± = q ± P/2, and Z x (by the Ward iden- 
tity Z\ = Z2) is the renormalization constant of the 



quark-photon vertex. The use of rainbow-ladder here en- 
sures that the important vector and axial-vector Ward- 
Takahashi identities (WTIs) hold [39, 40]. 

For our purposes it is reassuring to have a description 
that is able to reproduce a number of observables whilst 
at the same time being sufficiently simple that we can 
define unambiguously the hadronic four-point function. 
Furthermore, this approach describes at the same time 
both the perturbative and non-perturbative regime. This 
unification of IR and UV scales is very important for the 
problem of LBL and is not shared by effective low-energy 
descriptions of the strong interaction. The importance of 
this feature will be elaborated in section IV. 



III. CONTENT OF THE HADRONIC 
FOUR-POINT FUNCTION. 

Now we turn to the structure of the hadronic photon 
four-point function central to the computation 

of the hadronic LBL contributions. Using the rainbow- 
ladder truncation as introduced above, this object can be 
written exactly in the form (see Appendix B) 



(11) 



where the factors 6 and 12 indicate the number of per- 
mutations of the diagrams that one must consider. We 
see that there are two classes of diagrams: 

• The second class of diagram contains the T-matrix 
that describes all kinds of quark-antiquark inter- 
actions including the dynamical propagation of 
mesons. In our earlier work, the T-matrix was ap- 
proximated by pseudoscalar meson exchange [19, 
20] in reasonable agreement with low energy effec- 
tive models [11, 13-15, 17, 22, 41-46]. 

• The first diagram, which we refer to as the quark- 
loop topology, will constitute the main focus of this 
work. This object is composed entirely of fully 
dressed quark propagators and quark-photon ver- 
tices. In our previous publications [19, 20] we were 
not in a position to employ here the full quark- 
photon vertex as described by Eq. (10) due to the 
numerical complexity. Instead, we were limited to 
the Ball-Chiu construction [47] that fixes the first 
four components of Eq. (Al) exactly in terms of the 
quark dressing functions, see Eq. (A3). 

In this work we will investigate the leading transverse 
structure that, amongst other things, dynamically yield 
the picture of vector meson dominance (VMD) [32]. This 
is the case since these structures couple to the vector 
meson channel and hence one finds time-like poles corre- 
sponding to bound states. The leading component will 
be extracted and compared with its ENJL equivalent. 
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Note that we will limit the considerations to the case of 
two degenerate flavors for simplicity. The contributions 
of strange and charm quarks are only included for our 
best estimate of the LBL contribution in section VD. 

Note that there are contributions in the four-point 
function that are not accounted for in the representation 
of Eq. (11). These include unquenching effects due to 
internal quark lines that can be connected to the dynam- 
ical back-coupling of hadronic degrees of freedom [48] . In 
effective descriptions such contributions show up as pion- 
loops that organise themselves into a counting scheme 
within chiral perturbation theory. These contributions 
are typically considered to be sub-leading. A recent 
investigation, however, finds that next-to-leading order 
contributions might be more important than expected 
due to the somewhat accidental smallness of the leading 
terms [49]. 



IV. QUARK LOOP: COMPARING DSE TO 
ENJL 



In principle we have all of the ingredients at our dis- 
posal, see Eqs. (1,2) to calculate the quark-loop contribu- 
tion to LBL shown in Fig. 2. We note that this involves 
a three-loop integration over fully dressed quark propa- 
gators and quark-photon vertices which poses a consid- 
erable numerical challenge. Thus, to err on the side of 
caution we present a step- wise investigation of going be- 
yond the tree-level approximation reported in [19]. 

Additionally, it is useful to put our calculation in per- 
spective so that a greater understanding of our approach 
can be conveyed. To this end, we will present a compar- 
ison between the DSE approach and the, in many ways 
similar, ENJL model [50]. There are of course subtle 
differences between the two, and we will here point out, 
contrast, and discuss the consequences of each. This we 
further cement by testing ENJL-inspired vertices within 
our approach to explicate their kinematical differences. 




A. The ENJL perspective 

The inverse quark propagator in the ENJL model is 
just 



(12) 



} ENJL\ 

with p a Euclidean momentum and M the constituent 
quark mass. Note that the wave function renormalisation 
is just unity and the mass function is independent of 
momentum. 

The quark-photon vertex is given here by a bubble sum 



(13) 



see Refs. [11, 51] for details. Owing to the simple struc- 
ture of the effective four-quark vertices, the contributions 
can be resummed as a simple geometric series of one-loop 
amplitudes J2n Bubble™ = 1/(1 — Bubble). As a result, 
the quark-photon vertex depends only on the photon mo- 
mentum Q 2 and not on the relative momentum of the 
quarks 1 . 

Thus we have 



-iENJL 



(Q) = 7m _ .9ENjLn MI/ (Q)7 l/ 



(14) 



where (/enjl is derived from standard NJL couplings. 
The bubble sum n Ml/ (Q) has the transverse structure 



iWQ 2 ) = (qV - qM n T (o 2 



(15) 



since a potential longitudinal piece n L (Q 2 ) vanishes for 
identical quark masses [51]. 

According to Ref. [51] the transverse part in the VMD 
limit has the form 



n T (Q 2 ) = 



(16) 



where the momentum dependence of the mass function 
My and decay constant fy are neglected, which is re- 
ported to be a rather good approximation to the mo- 
mentum dependent case. In the limit we consider here, 

VI 



i v M v = 1/^enjl and so 



^ENJL 



7m 



T 

7m 



Q 2 



My 



(17) 



where 7J = 7 l /T Alt/ (Q) and My is identified with the 
mass of the p-meson. Together with the quark prop- 
agator Eq. (12) the vertex satisfies the vector Ward- 
Takahashi identity (WTI), Eq. (A2). Note that the ver- 
tex in Eq. (17) can be written as ~ 1^M V /{Q 2 + M v ) if 
the quark loop is transverse. This explicates why a sup- 
pression is found for VMD inspired transverse dressings 
compared to a bare vertex 7^. 



FIG. 2. The quark loop contribution to the muon g — 2. The 
quarks and vertices are dressed according to Eqs. (4,10). 



1 Note that a similar simplification happens for the four-quark 
T-matrix, which essentially reduces to a bubble sum sandwiched 
between two pairs of quark legs. The resulting mesons are point- 
like objects. 
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B. The DSE perspective 

In the DSE approach, the inverse quark propagator has 
the following covariant decomposition 



J DSEl 



(18) 



where in contrast to Eq. (12) we have a momentum de- 
pendent wave function Zf(p 2 ) and mass function M(p 2 ). 

In the rainbow-ladder truncation employed here, the 
quark-photon vertex is a dressed ladder-resummation of 
effective gluons 



I a 
I 
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$ 

C 






1 


E 


1 





(19) 



If we replace the gluon exchange by a momentum inde- 
pendent contact interaction we arrive at a picture similar 
to the ENJL model above. Keeping the exchange as is, 
the calculation is somewhat more involved as we can no 
longer have a simple geometric series. 
The full quark-photon vertex is given by 



fit, we will also employ the component of the vertex 
as extracted from the full calculation of the quark-photon 
vertex BSE. 



C. Differences between DSE and ENJL 

One of the key differences between the ENJL and DSE 
approaches is that in the former we have a contact in- 
teraction, whilst the latter features an interaction that 
features momentum exchange. This has far-reaching con- 
sequences as we discuss below. 

Firstly we look at the differences in the quark propa- 
gator, see Fig. 3. In the ENJL model we have Zf(p 2 ) — 1 
and M{p 2 ) = M const as opposed to the fully momentum 
dependent functions from the DSE. For the quark mass 
function, we see that in the DSE it saturates in the IR 
at about M(0) ps 450 MeV and continuously connects to 
its perturbative running at large momenta. The ENJL 
model, in contrast features a constituent-like quark mass 
of ~ 300 MeV at all scales up until the model cut-off 



ry q, *) = J2 A(i) 4 i) (Q. fc ) + E Tii)T P w> fc ) . ( 2 °) 



i=l 



i=l 



withAW = \^{Q 2 ,k 2 ,Q-k) and r« = r^{Q 2 ,k 2 ,Q-k) 
the longitudinal and transverse scalar coefficients respec- 
tively that correspond to the basis elements given in 

Eq. (Al). Additionally = Y^i X ® L ®(Q> k ) de " 
fines the Ball-Chiu vertex with the coefficients fixed by 
the vector WTI [47]. Our main interest here lies in the 

eight transverse components Xjl (Q,k) which couple to 
vector bound-states. For the sake of comparison with 
the ENJL model, we will take only the leading trans- 
verse component T^(Q, k) = 7J under consideration in 
this work. 

While we are mainly working with the full numerical 
result for the quark-photon vertex it is sometimes use- 
ful to also have a simple analytical form at hand, which 
captures the main features of the numerical solution. An 
approximate form for the leading component of the trans- 
verse part has been given in Ref. [32] which depends on 
both the relative and total momenta of the vertex 



r p (Q,fc)~r 



BC 



T w 4 AV 

— - 



h 



Q 2 



+ k 4 M v Q 2 



M 2 



(21) 



This form has been fitted to the full numerical solution 
for the quark-photon vertex obtained from its BSE. Here, 
lo and a describe the suppression of the amplitude for 
large relative and total momentum respectively. We find 
reasonable agreement with the numerical solution with 
lo = 0.66 GeV, a = 0.15 and N v f v /M v = 0.152, see 
Fig. 5. Additionally, Ny is a normalization factor. We 
consider this formula to be the VMD limit of the DSE 
quark-photon vertex. In addition to this BSE inspired 
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FIG. 3. Comparison of the mass function M and the wave 
function Zf for DSE and ENJL quarks. 
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~ 1 GeV. For the quark wave function renormalisation, 
we see that for a large momentum range Zf(p 2 ) < 1 
which constitutes a suppression of the quark propaga- 
tor with respect to the constant Z of the ENJL model. 
This can have several consequences for the quark-loop 
contribution to g — 2, Fig. 2, where there are four quark 
propagators. With Zf(p 2 ) < 1 in the DSE approach, 
we may expect a suppression of the contribution by a 
factor Zf(s) for some representative, 'average' momen- 
tum scale s. On the other hand, the momentum depen- 
dent quark-mass function allows for quark masses smaller 
than M(p 2 — 0) to be probed. Naively, this leads to an 
enhancement of the contribution. This may be equiva- 
lent to using a momentum-independent quark mass M(s) 
where s is some representative momentum scale. This 'ef- 
fective' mass may be surprisingly small and, indeed, such 
small quark masses have been observed to be necessary 
in several models [15, 17, 52]. We come back to this point 
below. 

If we now compare the ENJL vertex, Eq. (17) and its 
DSE equivalent, Eq. (21) we see that there is a similar 
structure. That is, there is a part dictated by the WTI 
(gauge part) and a part that represents the VMD physics 
of the transverse vertex. Whilst in the case of the DSE we 
have functions that depend on both the photon and rel- 
ative quark momenta, in the case of the ENJL model we 
have trivial momentum dependence for the gauge-part, 
and reduced momentum dependence for the transverse 
part. 

First, let us discuss the gauge part of the quark-photon 
vertices. The leading coefficient of 7 M has the form 



X^(Q 2 ,k 2 ,k-Q) 



+ 



1 



Z } (k\) Z f (k 2 _) 



(22) 



where k± — k ± Q/2. Based on the difference in be- 
havior of the quark propagators, we see for the ENJL 
model that Ai = 1, whilst for the DSE Ai > 1. Thus, 
for the DSE we expect the gauge part of the vertex to 
yield an enhancement of the quark- loop ~ 1/Zf(s) 4 as 
to the bare vertex approximation. A comparison of these 
components of the vertex is shown in Fig. 4. 

Now we take a look at the dominant transverse compo- 
nent of the quark-photon vertex. In Fig. 5 we show the 
dressing function as calculated sclf-consistcntly from the 
DSE for the quark-photon vertex, together with the fit 
function of Eq. (21) and the equivalent part of the ENJL 
vertex, Eq. (17). We see that the fit to the DSE, as 
a function of the total momentem Q, behaves very simi- 
larly as the full numerical solution in the dominant region 
around the scale My = 0.77 GeV. Small deviations occur 
at large momenta due to the exponential fall-off of the fit 
function. The dressing function of the ENJL model has a 
similar behavior except at large momenta where no fall- 
off is seen and instead it tends to a constant. However, 
due to the weighting of the integrand in the calculation 
of it transpires that such deviations at large Q 2 are 
not relevant. 



However, differences between the DSE and ENJL ap- 
proach become readily apparent when one considers the 
impact of the relative momentum on the kinematics. 
In the DSE approach the transverse VMD-piece is sup- 



Aq CD , an effect which is 



pressed for momenta k 2 > 

not present in the ENJL approach due to its contact in- 
teraction. As a consequence, we expect smaller negative 
contributions to LBL in the DSE approach as compared 
to ENJL. The degree of overestimation of negative contri- 
butions in ENJL depends upon the kinematical weighting 
of the integrand, and can be estimated in our approach 
by systematically removing the dependence on the rel- 
ative momentum in our quark-vertex dressing function. 
This will be detailed in the next section. 



V. RESULTS 

Here we present a more quantitative analysis of the 
difference between the ENJL and DSE approaches, with 
respect to various approximations that can be made in 
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FIG. 4. The leading A' 1 ' component of the quark-photon ver- 
tex constrained by the Ward-Takahashi identity of Eq. (A2) 
and given in Eq (22). We show two slices as a function of 
the total and relative momenta Q 2 (with k = 0) and k 2 (with 
Q = 0), respectively. Note that the constant dressing corre- 
sponds to the ENJL model. 
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FIG. 5. The dependence of the dominant transverse dressing 
on the photon momentum Q 2 is shown for the explicit 
solution of the quark-photon vertex BSE, Eq. (10), the fit to 
this given in Eq. (21), and the transverse part of the ENJL 
vertex, Eq. (17). 



the calculation of the quark loop. In order to compare 
the two approaches we need a means to find the most im- 
portant and therefore representative momentum scales. 
This we achieve by averaging over the dressing functions, 
as weighted by the importance sampling of the VEGAS 
Monte-Carlo [5-3] we use to evaluate the quark-loop con- 
tribution to LBL. 

For the calculation in the DSE-framework we use an 
IR-cutoff of 1(T 3 GeV and an UV-cutoff of 10 3 GeV for 
all three loop variables. No additional splitting into a 
perturbative and a nonperturbative region is necessary, 
since these regions are continuously connected, see the 
discussion around Fig. 3. In contrast, this is different for 
the ENJL model where the calculation is split into a low 
and a high-energy part [1 1]. In the former, constituent 
quarks Eq. (12) with M = 0.3 GeV are used together 
with the vertex construction Eq. (17). This low-energy 
part is cut off by requiring the photon momenta to be g < 
H, where /i is varied in the range between 0.7 GeV and 4 
GeV. The high-energy contribution is approximated by 
a bare quark-loop (bare propagators and bare vertices) 
where the mass of the quarks is given by the same scale 
H which is supposed to act as an effective IR-cutoff. It 
was found in [11] that the result is rather stable against 
a variation of fi in the considered range. The final result 
is quoted as 



a 



LBL, ql, ENJL 



21(3) x 10 



-li 



(23) 



Note that if fj, = 4 GeV is taken, almost the complete 
contribution comes from the low-energy part. While this 
value is usually considered to be inappropriately high for 
an NJL type model, we find this viewpoint useful since it 
facilitates the comparison to the DSE case, where there 



Quark dressing 


(Zf) 


V(M 2 } [GeV] 


Mio- 11 ] 


Zf and M dynamical 


0.83 


0.14 


47 


Zf = Z% , M dynamical 




0.18 


102 


Z t = Z 2 -\ M = 0.477 






22 


Zf = Zf 1 , M = 0.300 






51 


Zf = Z2 1 , M = 0.200 






104 



TABLE II. The quark-loop contribution to hadronic LBL to- 
gether with effective average values for the Zf and M func- 
tions determined with a bare quark-photon vertex /?27 M . We 
compare results for dynamical quark dressing functions to 
those with static ENJL-like equivalents. The quark wave- 
function renormalisation constant Z2 ~ 0.97 for the cutoff 
and renormalisation point used. 



is no separation in high and low-energy parts. We will 
later see that the contribution Eq. (23) can be roughly 
reproduced even with the much higher cutoff we use in 
our calculations 2 . The results presented in this work are 
obtained within Landau gauge for QED (£ = in the 
photon propagator apparent in fig. 2), in contrast to our 
earlier work [19, 20] where we used Feynman gauge for 
QED (£ = 1). It turns out that all the results presented 
here are affected only mildly by this choice, as can be 
scon in Table VI. 



A. Impact of a dynamical Zf and M 

To gauge the impact of momentum dependent quark 
dressing functions, we compare several calculations in 
the bare-vertex approximation, see Table II. Note that 
we attach a renormalisation constant to the bare vertex 
e.g. Z 2 'y l _ l to ensure multiplicative renormalisability. For 
a momentum independent quark wave-function we also 
supplement the quark by a factor Z^ 1 for the same rea- 
son. These extra factors cancel in all but the first line of 
Table II. The corresponding case of bare vertex with fully 
dressed quark thus differs from our earlier results [19, 20] 
by a factor of Z% ~ 0.89 in addition to the different QED 
gauge and the smaller number of quark flavours. 

Let us first compare the first two lines of Table II. Re- 
placing the dynamical quark wave function by Zf = Z^ 1 
we note an enhancement of the contribution to a M roughly 
by a factor of 2. This is in good agreement with the 
expectation discussed above: since there are four quark 
propagators in the quark loop, the full calculation con- 
tains an extra factor of the order (Zf) 4 with (Zf) the 
average value of the wave-function that is probed. In- 



deed, one finds (Z f y ~ (0.83) 4 - 0.5. 



Note furthermore that Eq. (23) includes strange quarks which 
we neglect in this work. These give, however, a few percent 
correction [1 1] such that this detail is not important for the more 
general analysis we present here. 
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Quark Dressing 


(Zf) 


a^lO" 11 ] 


Z/ and M dynamical 


0.76 


100 



TABLE III. The quark-loop contribution to hadronic LBL 
with the 1BC dressing and the average of the Zf function. 
These are to be compared with the bare-vertex results in the 
previous table. 

To explore the impact of the momentum dependent 
quark mass function, we keep Zf = Z^ 1 fixed and com- 
pare several constant values of M to the dynamical case. 
For M — 0.477 GeV, the infrared plateau of the dynam- 
ical mass function, we obtain less than l/5th of the dy- 
namical result (with Zf = Z^ 1 )- A quark mass of M = 
0.3 GeV gives a result commensurate with the calculation 
of Ref. [ ] which yields a ^ BL >i l > N f= 2 = 49.1(3.4) x lO" 11 
for the contribution of bare u and d quarks of the same 
mass M = 0.3 GeV 3 . Only a significantly smaller con- 
stant mass ~ 0.2 GeV close to the average one probed un- 
der the integrand leads to a result that compares with the 
dynamical one. Note, however, that this does not work 
in general, i.e. the fully dynamical result of 47 x 10~ n is 
not reproduced by a static approach with Z — 0.83 and 
M = 0.14 GeV. 

B. Impact of dressed vertices: Gauge Part 

Here we focus upon the gauge part of the vertex. In 
the case of the ENJL model this part is just bare since 
Zf = 1. In the DSE approach for the leading part, 
Eq. (22), we give the fully dynamical result in Table III. 
By comparison with the first value of Table II we find 
an increase by more than a factor of 2 when the ver- 
tex dressing is included. Again, this is roughly what one 
expects, since the contribution from the vertices gives a 
factor ~ (Zf)~ 4 . Note that the enhancement from the 
gauge-part of the vertex is comparable to the suppression 
due to a non-trivial Zf in the quark propagator. 

In principle, we should include not only the leading 
term of the gauge part of the vertex but the full Ball-Chiu 
vertex as fixed by the Ward-Takahashi identity. However, 
this is currently not possible due to numerical instabilities 
which are currently not under control 4 . 

C. Impact of dressed vertices: VMD Part 

To better compare with the ENJL model we consider 
first the case where the quark-photon vertex is taken to 
be of the form 

r M ( Q, *) = Z 2l , + r«(Q, k) T« (Q, k) , (24) 



3 This number is extracted from Table I of Ref. [54] 

4 The corresponding results given in Ref. [20] are presumably not 
correct as will be detailed in an erratum. 



with the leading transverse component T^\Q,k) = 7J, 
and T^(Q,k) its dressing function. The quark wave 
function renormalisation constant Z 2 is necessary to 
maintain multiplicative renormalizability. For the trans- 
verse dressing function t^'(Q, k) we study three choices: 
(i) the ENJL model, Eq. (17), which we also supplement 
by an additional factor 

f(k 2 ) = l/(l + fc 4 /0.66 4 ) (25) 

to simulate a relative momentum dependence, (ii) the 
DSE approach with the numerical solution to Eq. (10) 
and (iii) the DSE approach with the analytic fit func- 
tion Eq. (21) to the numerical result. Both, (ii) and (iii) 
already include a relative momentum dependence in con- 
trast to the original ENJL approach. Our results are 
presented in Table IV. 

First of all note that the result where a bare vertex 
is used together with the ENJL type transverse part 
(first line) is quite close to the ENJL result shown in Eq. 
(23). We are thus able to reproduce the ENJL result 
numerically, despite the very different UV cutoff, that 
obviously does not matter much. Taking a closer look at 
Table IV we see a common pattern. When the vertices do 
not depend on the relative momentum between the two 
quarks, the contribution is around ~ (14—16) x 10 -11 . 
Compared to the corresponding result with bare-vertex 
of 47 x 10~ n we therefore find a reduction of similar size 
as the one from ~ 60 x 10~ n to ~ 20 x 10 -11 reported 
by [11]. We obtain this large suppression in our ENJL 
type calculation as well as in the DSE approach with mo- 
mentum dependent quark propagators. However, when 
we take into account the relative momentum dependence 
in the full DSE calculation and including the additional 
factor f(k) in the ENJL model we find that the suppres- 
sion due to transverse parts is much reduced. We find 
results in the range of a M ~ (41 — 46) x 10~ . This is 
at most a reduction of just ~ 15% and constitutes one of 
the main results of this work. 

Finally, we give the same comparison with the inclu- 
sion of the leading L\ dressing \^ from Eq. (22) given 
by the WTI, see Table V. We see that the trend here 
is very similar as for the bare vertex, except now with 
an enhancement due to the non-trivial dressing function 
of the gauge part. This enhancement is of the order of 
2 — 3. It differs slightly for the cases with and without 
relative momentum dependence, which shows that there 
is also interference between the different vertex compo- 
nents. Note that the A^ 1 ) dressing is always used with its 
full kinematics. 



D. Best result and electromagnetic gauge 
invariance 

Our most reliable estimate for the contribution from 
the quark-loop to the anomalous magnetic moment of 
the muon is the one obtained with full dynamics in the 
quark propagator, the leading gauge-part A*- 1 -* and the 
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Vertex Dressing 


a^lO -11 ] 




14 


n T r 1 i 2 \ 

^27m +7m t enjl/(k ) 


45 


^2 7 M+7j^'(fc = 0) 


16 


Z 2 7 M + 7 Jr« 


46 


^7M+7jri 1 1 ' c (fe = 0) 


14 


^27/, + 7^1^ 


41 



TABLE IV. Bare and leading transverse vertex component, 
with dressing functions from the ENJL model, VMD like fit 
from DSE/BSE (t^), and from an explicit calculation of the 
quark-photon DSE (t£j^). Results are shown with, and with- 
out the inclusion of a dependence on the relative momentum. 



Vertex Dressing 


a M [10- n ] 


7m^ (1) + 7j r ENJL 


43 


7^ (1) +7jr- ENJL /(fc 2 ) 


103 


7,A«+7j^'(fc = 0) 


43 


7„A«+7?r« 


105 


7mA (1) + 7jr c «(fc = 0) 


41 




96 



TABLE V. 1BC and leading transverse vertex component, 
with dressing functions from the ENJL model, VMD like fit 
from DSE/BSE ("TgJ ), and from an explicit calculation of the 
quark-photon DSE {t^{ c ). Results are shown with, and with- 
out the inclusion of a dependence on the relative momentum. 



Interaction 




a LBL,QL x 10 11 


MT 







96 


MT 


l 




94 



TABLE VI. Results for the lBC+transverse vertex dressings. 
We compare different photon gauge parameters £ = {0, 1}, 
i.e. Landau and Feynman gauge. 

\(i..3) are necessary to maintain electromagnetic gauge 
invariance. As mentioned above, this is currently not 
possible due to severe numerical problems with the terms 
A^ 2 - 1 and A^ 3 ). We therefore have to gauge the error in the 
present calculation due to violations of gauge invariance. 
This is conveniently done by varying the QED gauge 
parameter £. Results for Feynman and Landau gauge 
are shown in table VI. The variations with £ are on the 
two percent level and therefore reassuringly small. The 
insensitivity with respect to the gauge parameter is an 
indicator for the (almost) transversality of the resulting 
quark- loop part of the photon four-point function 11^,^. 

We emphasize, however, that the smallness of the 
gauge violations due to the omission of A' 2 ) and A' 3 ) 
cannot be taken as an indication that these terms will 
not contribute much to the physical, transverse part of 
the photon four-point functions. As already mentioned 
above, our previous calculation of these contributions in 
Ref. [20] is presumably not correct and needs to be thor- 
oughly reinvestigated. This will be done in future work. 



leading transverse part t^] of the quark-photon vertex. 
For two light quark flavors we obtained 

LBL,Ll+Tl,iVf=2 /n „ . „s 1n -ll ,r, c \ 

' = (96±2)xl0 , (26) 

where the error is purely numerical. Compared to the 
corresponding value a^ BL ' L1 ' r ~ 2 = (100±2) x 10~ n for 
the 1BC vertex without transverse parts we thus find a 
suppression of the order of five percent due to the VMD 
physics. This is much less than in simple models. If we 
additionally include the strange and charm quark contri- 
butions we arrive at 

a LBL,Ll + Tl,A f/ =4 = ^ ± 2) ^ 1Q _u (2?) 

which compares to the 1BC Nt = 4 case (in Landau 
gauge) a L BL . L1 ^/= 4 = (in ± 2) x lO" 11 . Note once 
more that the corresponding result in [21 1] differs slightly 
due to QED Feynman-gauge, see discussion below. 

We are currently working on the further inclusion of 
the other transverse terms r' 2 " 8 \ corresponding results 
will be presented elsewhere. The potential impact of 
these terms is hard to gauge without an explicit calcula- 
tion. Nevertheless, from a systematic point of view, the 
omission of these terms is unproblematic. This is differ- 
ent for the non-transverse part of the vertex. Strictly 
speaking, the presence of all three Ball-Chiu components 



VI. SUMMARY AND CONCLUSION 

In this work we investigated in detail key similarities 
and differences between the ENJL and DSE approaches. 
Whereas the ENJL model features a contact interaction 
giving rise to a trivial quark wave-function renormalisa- 
tion and a constant quark mass function, in the DSE ap- 
proach these are both momentum dependent. The same 
is true for the gauge part of the vertex, that is determined 
by a Ward-Takahashi Identity. Whereas in the ENJL 
model this part is trivial, the corresponding Ball-Chiu 
terms in the DSE approach are nontrivial and momen- 
tum dependent. In both approaches there are transverse 
parts in the vertex which are dominated by the vector 
meson poles leading to a characteristic behavior also in 
the space-like momentum region. 

When assessing the influence of the different momen- 
tum dependent dressing functions of the DSE approach 
as compared to ENJL we found partial cancellations. 
Dressing effects due to the non-trivial wave-function on 
the level of 50% are cancelled by opposite effects due to 
the dressing of the gauge part of the vertex. An im- 
portant effect that is not cancelled is the one of the dy- 
namical mass function. We found that this function is 
not tested predominantly at its large infrared plateau of 
M(0) ~ 0.48 GeV but rather at smaller values at inter- 



10 



mediate momenta commensurate with M ^0.2 GeV. 
This gives rise to a larger contribution than expected 
from constituent quark-loop calculations. Our finding 
may serve to explain the surprisingly small (constituent) 
quark masses needed in chiral models to obtain sensi- 
ble results [17]. Our most interesting result, however, 
is related to the transverse part of the vertex. These 
are of high interest because they dynamically include the 
phenomenology of vector meson dominance and are ex- 
pected to reduce the overall light-by-light contribution. 
As in the ENJL model the corresponding behavior of the 
transverse part of the vertex in the DSE-approach is gen- 
crated dynamically. A key difference is, however, that in 
the DSE case the dependence on the relative momen- 
tum of the quarks is taken into account. On the level of 
mesons (as e.g. the p-meson that is the vital ingredient 
in VMD) this takes the distribution of quarks inside the 
bound state into account, which is not the case in the 
ENJL model [55]. The inclusion of this effect, either in 
the full DSE calculation or by a suitable modification of 
the VMD-term in the ENJL model has important conse- 
quences: the reduction of the quark-loop contribution to 
due to VMD effects, observed in previous calculations 
[ ], is drastically reduced. 

Thus we can pinpoint the differences between the DSE 
and ENJL calculations to be down to the contact inter- 
action limiting the momentum dependence of dressing 
functions. We believe that ignoring the relative quark 
momentum in the quark-photon vertex overestimates the 
suppression that the transverse part of the vertex pro- 
vides and thus lowers significantly it's numerical contri- 
bution to the anomalous magnetic moment of the muon. 

Our present best result for the quark-loop contribution 
has been discussed in section VD. Combined with our 
result for the pseudoscalar meson exchange diagram from 



Ref. [20], 
estimate 



-,LBL:PS 



LBL 



(81 ± 2) x 10~ n , we arrive at the 



(188 ± 4) x 10" 



(28) 



for the total LBL contribution. Again, the error is purely 
statistical due to our numerics. Since at present any 
guess of the systematic error of this number is clearly 
subjective (due to the omission of terms in the quark- 
photon vertex), we do not attempt such an estimate. 

It has to be emphasized, however, that our determi- 
nation of the quark loop contribution to light-by-light, 
Eq. (28) is by far not complete since several terms in the 
vertex dressing are still missing. The study if the influ- 
ence of these terms is an important task for the future. 
Nevertheless we hope that the systematics of the present 
work serves to give the reader a better understanding of 
the technical and physical mechanisms at work in the 
complicated light-by-light scattering contribution. Fur- 
thermore we have shown that our best results at present 
are stable under variations of the photon gauge parame- 
ter and therefore serve as an important intermediary step 



towards a full calculation to come. 
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Appendix A: Quark-photon vertex 

In this appendix we give the explicit basis for the 
quark-photon vertex. It was taken from Ref. [ ] . 



L\ = 2k^ 



-3 

L 4 



K = 2ifc M 



T 1 = 7 2 



T 

H — in 



1 



Tl^[klx T -^il{k T fye 
it = k l? p ■ k /( k2p2 ) 

Tt = -( 1 T\P ) X T } + 2k^P)/2k 
(P-kf 



(Al) 



t; = ihl,A ( i 



" v 1 V P 2 k 2 

TisT ml 1,2 



2T° 



Here P is the total momentum (photon momentum) and 
the relative momentum is k = (k + + /c_)/2 were k± are 
the quark momenta. The symbol T denotes transversal- 
ity with respect to P. Furthermore, k = \Jk 2 . The twelve 
dressing functions corresponding to the tensor structure 
(Al) are denoted as A^, (i — 1, 2,3,4) corresponding to 
the and with (j — 1, ... ,8) corresponding to the 

T/f', see Eq. (20). The dressings \^ are determined by 
the WTI 



iP^P, k) = S- 1 (k.)- S-\k + ). 



(A2) 



and regularity demands [ ]. The resulting vertex con- 
struction is referred to as the Ball-Chiu construction 



rjf (P, k) = J^A + tyk^A + ifc^As, 



(A3) 



where the symbols 

F(k 2 + ) + F(k 2 _) 



_ F{kl)-F{kl) 

A F - f,2 _ k 2 . ( A4 ) 



have been used and A and B are the quark dressings. 
Note that A' 4 ) is identically zero. 
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FIG. 6. The photon DSE. Full Propagators are denoted by 
white blobs and the red blob marks the full 1PI fermion- 
photon vertex. 



Appendix B: Derivation of the hadronic four-point 
function 



In the present appendix we present the derivation of 
the photon four-point function that is needed in the light- 
by- light scattering contribution to the muon g — 2. The 
four-point function is obtained by taking two derivatives 
of the full inverse photon propagator 



n 



S 4 T[A] 



SA^AvSAaSAp 

5 2 S 2 T[A] 5 2 
SA^SAv SA a 6A p ~ SApSAv 



(Bl) 



In order to proceed we consider the photon DSE that is 
presented in Fig. 6. We now take two further derivatives 
of the photon DSE. By the virtue of the approximation 
we consider, namely rainbow-ladder truncation of QCD, 
the photon is not a dynamical field in our formalism. 
The photon self-energy that we calculate is performed 
using quark and quark-photon vertex equations that do 
not take into account the photon as a dynamical object, 
but rather as an external background field. Internally 
all objects just include dynamical quarks and model 
gluons. This is why the application of further derivatives 
on the photon DSE is consistent within our approxi- 
mation. Furthermore this operation leaves the internal 
consistency of the truncation intact. In particular the 
consistency with chiral symmetry and electro-magnetic 
current-conservation in form of WTI's is not destroyed. 
This technique is also referred to as 'gauging' [5f ]. For 
a detailed description of this procedure in particular in 
rainbow-ladder truncation see [57]. 

An object that will be needed on several occasions 
is the quark anti-quark photon-photon vertex. Within 
rainbow-ladder approximation this object can be exactly 
written in terms of quark-photon vertices and the T- 
matrix [57] 




= 2 



(B2) 



where the factor of two is not a symmetry factor but 
rather a convenient way to denote that the two possible 
photon permutations are included. Application of one 



derivative to the photon DSE in Fig. 6 gives 






The bare inverse photon does not contribute. The deriva- 
tive has to be applied to the quark-photon vertex and the 
two dressed quarks of the photon self energy. Since there 
are two quarks in the loop the second diagram comes 
again with two permutations that correspond to one dia- 
gram with straight photon legs and one with crossed ones. 
This is signaled by the factor 2. Now the application of 
the relation shown in Eq. (B2) results in 





T 



Next we take into account the relation between the T- 
matrix and the G-matrix 



G — SS -\- S ST S S , 



(B3) 



which says nothing but that the T-matrix is the ampu- 
tated connected part of the G-matrix. In Eq. (B3), S is 
the quark propagator and multiplication indicates con- 
traction of dirac-indices and integration over momentum 
arguments. The inhomogeneous part that involves just 
two propagators (SS in Eq. (B3)) cancels the second di- 
agram on the right side of the equation above and we are 
left with 






Now the final step is to remember that the bare quark- 
photon vertex together with the G-matrix corresponds to 
a dressed quark-photon vertex and to external quark-leg 
dressings. Finally we obtain the consistent representa- 
tion of the hadronic part of the three-photon vertex 



= 2 



which is not zero as long as the external photon 
background- field is present. The factor 2 can also be 
interpreted as representing the two possible orientations 
of the quark loop. 

Now the four-point function is obtained by taking a 
further derivative with respect to the external photon 
field. Again the derivative can act on quark-photon ver- 
tices and the quark propagators. The result is 
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where the additional factor three is caused by the three 
possible vertices and three propagators respectively. Us- 
ing again relation (B2) we arrive at 




Up to this point no approximations on top of rainbow- 
ladder truncation of QCD have been made. The last two 
representations of the four-point function are thus truly 
consistent with the rainbow-ladder photon self-energy 
shown in Fig. 6. 
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